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ABSTRACT 



We show that the effective coarse graining of a two-mode squeezed density 
matrix, implicit in the Wehrl approaches to a semiclassical phase-space 
distribution, leads to results in agreement with previous different defi- 
nitions of entropy for the process of pair production from the vacuum. 
We also present, in this context, a possible interpretation of the entropy 
growth as an amplification (due to the squeezing) of our lack of knowledge 
about the initial conditions, which gives rise to an effective decoherence 
of the squeezed density matrix. 



Some years ago it was suggested that the in quantum process of pair production 
from the vacuum, induced by the action of a time-varying classical background, 
the growth of the average number of quanta of a given field should be naturally 
associated to the growth of entropy of that field [1,2]. More recently, such an 
entropy growth was quantified on the grounds of a squeezed state approach to the 
process of particle production [3-7]. With this approach one can easily compute, in 
particular, the total entropy S stored in the cosmological perturbations amplified 
by inflation, and one obtains [4,6] 



where S C mb is the thermal black-body entropy of the 3K° photon microwave back- 
ground, Hf the final curvature scale of inflation, and M p the Planck mass. 

The use of the squeezed state formalism to compute entropy production might 
seem to give rise to a puzzle, since the squeezed evolution is unitary. The solu- 
tion of the puzzle is that, of course, there is no loss of information in principle 
associated to the evolution of the initial state into a final squeezed quantum state 




(1) 
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(and so no possible information paradox, like in the case of black-hole dynamics). 
The loss of information occurs however in practice, as a consequence of the process 
of measure of the observables characterizing the final squeezed state. There are 
indeed variables (related to the so-called "superfluctuant" operators [8]) which are 
amplified by the squeezed evolution and become thus available for macroscopic ob- 
servations, while other variables (the "subfluctuant" ones [8]) are correspondingly 
"squeezed" and become unobservable. 

A simple and intuitive way to illustrate this situation is to consider a semi- 
classical description of the squeezing process as a parametric amplification of the 
wave function, "hitting" the effective potential barrier of a Schroedinger-like equa- 
tion (with time-like variable) [9,10]. The final amplitude of the wave scattered by 
the barrier may be larger or smaller than the initial amplitude, depending on the 
phase with which the wave enters the barrier. What is macroscopically detected, 
however, is the average amplification of the wave, out of an initial random distri- 
bution of phases. This point was clearly stressed already in the first paper on the 
cosmo logical amplification of tensor perturbations [11]. A macroscopic observa- 
tion of the quantum fluctuations of the metric, amplified by the cosmic evolution, 
traces out over phases. 

In order to compute the entropy associated to a squeezing process, this macro- 
scopic loss of information has to be taken into account in the form of a suitable 
"coarse graining" of the density matrix of the final (pure) squeezed state. Such a 
coarse graining was originally performed in a Fock basis, i.e. with respect to the 
eigenstates of the number operator, by averaging over phases [3-5]. The number 
operator, however, is only a particular example of superfluctuant variable which 
grows in the case of a pair production process, while the conjugate variable (the 
phase) is correspondingly squeezed. This suggests the use, in general, of the eigen- 
states of the superfluctuant operator as the natural basis for the coarse graining of 
a squeezed density matrix [6,7]. A possible different basis, that of coherent states 
with real amplitude, has been suggested however in [12]. 

It is important to stress that all these different approaches [3-7,12] to the 
entropy of a squeezed state exactly lead to the same result when the entropy 
growth ASk, for each mode k, is evaluated in the large squeezing regime, namely 

AS k ~ 2r k , r k » 1 (2) 

where r k is the modulus of the (complex) squeezing parameter [8]. In the previ- 
ously quoted papers, however, the decoherence process which leads to a reduced 
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density matrix p re d (from which to compute the quantum Von Neumann entropy) 
is imposed "by hand" , by setting the off-diagonal terms of the matrix to zero in 
some chosen basis; at best, like in the case of the number basis, it is justified in a 
semiclassical way using a stochastic approach and a sort of random-phase approx- 
imation [3,4]. The problem is that it is formally impossible in the number basis, 
as well as in the more general superfluctuant basis, to compute the quantum von 
Neumann entropy (—Trp re d In p re d) from a reduced matrix (p re d ^ P^ed) obtained 
by tracing out phases, or another subfluctuant variable. Indeed, a subfluctuant 
operator does not commute with the conjugate superfluctuant one, so that both 
operators cannot be simultaneously included in the complete set of observables 
characterizing a pure squeezed state. 

In view of these preliminary remarks, we can say that the aim of this paper is 
twofold. First of all we want to compute the semiclassical Shannon- Wehrl entropy 
[13,14] for a two- mode squeezed vacuum state, in order to show that also in that 
approach one recovers the result (2) for the large squeezing limit. This confirms 
that the entropy of a highly squeezed state is very robust to the particular scheme 
of coarse graining implemented (as discussed also in [15,16]). 

In the second place we want to present a possible justification to the reduction 
of the density matrix in a superfluctuant basis [6,7], by interpreting such a reduc- 
tion as the consequence of our lack of knowledge about the initial configuration 
which is the starting point of the squeezed evolution. To this purpose we shall 
represent the initial state, in general, as a statistical mixture of coherent states. 
This representation automatically defines a quantum von Neumann entropy asso- 
ciated to the initial mixture. This entropy however is constant, as the squeezed 
evolution is unitary, and then it does not contribute to the difference AS between 
final and initial entropy. However, if the statistical distribution satisfies (at least 
approximately) a principle of equipartition of the probabilities among the states 
of the mixture, then the off-diagonal terms of the density matrix turn out to be 
suppressed. An effective decoherence is thus produced, whose associated entropy 
in the superfluctuant basis is eventually amplified by the squeezed evolution. 

The arguments presented in this paper are intended to apply to any process of 
pair production, described by a Bogoliubov transformation connecting the particle 
{6, tf} and anti-particle {&, b'} annihilation and creation \in) operators to the \out) 
ones, {a, a\ a, a^} . Such a transformation reads, for each mode k, 

a k = c + (k)b k + c*_(k)P_ k 

al k = c_(k)b k + c* + (k)bl k (3) 
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where the Bogoliubov coefficients c±(k) satisfy |c + | 2 — |c_| 2 = 1. By putting 

c+(k) = coshr fc , c*_(k) = e 2i0k sinhr fc (4) 

the transformation (3) can be rewritten as a unitary transformation, 

a k = T, k b k T,\, a j _ k = Y< k V_ k Y\. (5) 

generated by the two-mode squeezing operator [8] 

E fc = exp (z* k b k b_ k - z k b{b^ , z k = r k e 2i0k (6) 

Our discussion can thus be applied in general to any dynamic situation in which the 
initial vacuum state |0) evolves towards a final two-mode squeezed vacuum state, 
| z) = Sfc|0), characterized by a non- vanishing value of the squeezing parameter 
r k = \z k \ (related to the expectation value of the number of produced particles by 
(n k ) = {z k \b\b k \z k ) = (0\a[a k \0) = |c_(/c)| 2 = simi 2 ^). 

It is important to note that in the context of such a dynamic evolution one can 
always define two operators x and x, called superfluctuant, whose variances (Ax) z , 
(Ax) z are amplified with respect to their vacuum value (Ai)o, (Ax)q, namely [8] 

(Ax) -> (Ax) z = (Ax) e r , (Ax) -> (Ax) z = (Ax) e r (7) 

where {Ax) 2 z = (z\x 2 \z) — ((z\x\z)) 2 (and the same for x). The variance of 
the canonically conjugate operators y, y is correspondingly squeezed, (Ay) z = 
(Ay) e~ r = (Ay) z . In terms of the superfluctuant variables the operators b and 
b^ have the differential representation [6,7,8] 

b = l -e ld (x -ix + d x - id s ) 

b= -e i0 (x + ix + d x +id s ) (8) 

(henceforth the mode index k is to be understood, if not explicitly written, and 
any correlation among modes with different \k\ will be neglected, following the 
coarse graining approach of [2]). The normalized wave function t(j z for the two- 
mode squeezed vacuum (such that aifj z = = aip z ) becomes, in the (x,x)-space 
representation [6,7], 

1/2 

ip z (x,x) = (xx|E|0) = (xx\z) = ^-j e ~%( x2 + x2 ) j 
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a = e" 2r , (z\z) = 1 (9) 

The squeezed evolution |0) — > \z) = E|0) may thus represented, in the superfluc- 
tuant basis, scaling transformation — > y^x, with real positive 
parameter cr < 1, related to the squeezing parameter r according to eq.(9). 

In order to compute the Shannon- Wehrl entropy for a squeezed state we recall 
that the classical entropy is defined, for any given phase-space distribution f(q,p), 
as 

Sd(f) = - J dqd P f(q,p)\nf(q,p) (10) 

while the quantum (or von Neumann) entropy is defined in terms of the density 
operator p as 

S q (p) = -Trplnp (11) 

For a pure state, in particular, p 2 = p and S q = 0. Suppose now that we want to 
compute, for a given quantum state represented by p, its entropy in the "semiclas- 
sical" limit. It is well known that a possible procedure to obtain the semiclassical 
limit of a quantum observable is to compute the expectation value of that observ- 
able in the basis of the coherent states \a). According to Wehrl [14], a natural 
semiclassical phase-space distribution associated to p is thus given by 

f(a) = -(a\p\a) (12) 

(the so-called Glauber distribution function [17]), which leads to define the semi- 
classical (Wehrl) entropy Sw as 

Sw = S cl (/(a)) (13) 

For a two-mode squeezed state \z) (see [18,19] for the one-mode case) we have, 
in particular, 

1 

f(a, a) = — ~{aa\z){z\aa) (14) 
where the coherent states \aa) satisfy 

b\aa) = a\aa), b\aa) = a\aa) (15) 

(the operators b, b are defined in eq.(8); their corresponding eigenvalues are com- 
plex numbers, a = a\ + i«2, a = ct\ + ia-i). By expanding in the superfluctuant 
basis \xx) one easily obtains 



(xx\aa) = —= exp 

'7T 



-(3f^) 2 - -(Kfc) 2 + kx + ~kx--(x 2 + x 2 ) 
2 2 2 
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(aa\aa) = 1 (16) 
where 5fc denotes the real part of z, and 

k = -ie~ ld (a + a), k = -e~ ie '(a - a) (17) 



Moreover 



f(a,a) = — g y dxdxdx'dx'ip z (xx)ipl(x'x')(aa\xx)(x'x' [ 



aa.) = 



= 1 exp [-tanhr(aae- 2ie + a*a*e 2 * e ) - |a| 2 - |a| 2 ] (18) 
7T 2 cosh r 

With this Glauber distribution, the semiclassical Wehrl entropy (13) for the two- 
mode squeezed vacuum is then 

Sw = — J d 2 ad 2 af(a, a) In f(a, a) = 2 + 2 In n + 2 In cosh r (19) 

The same result holds for a squeezed- coherent state \zaoao) = £|cto<So), obtained 
by applying a squeezing transformation to the "displaced" vacuum |ao«o) = 
D(ao)D(6to)\0), where D(a) = exp(a&t — a*b) is the Glauber displacement oper- 
ator. 

It is important to note that for r >> 1 this semiclassical entropy reproduces 
the large squeezing behavior of eq.(2), previously obtained [1,3-7,12] with different 
approaches to the entropy of a particle production process. We also note that the 
entropy (19) is exactly the sum of two "marginal" entropies, 

S w = S™ &1 + S^~ 2 (20) 

obtained by decomposing the probability distribution /(a, a) in the complex a, a 
planes, and by tracing out one of the two cartesian variables, 

S^a^-f da 1 da 1 f(a 1 ,a 1 )\nf(a 1 ,a 1 ), f(a 1 ,a 1 ) = J da 2 da 2 f(a,a) 

S ™a 2 = ~ J da 2 da2f(a2,a2)lnf(a2,a 2 ), f (012,612) = J da 1 da 1 f(a,a) (21) 

This is a non-trivial result, as it implies that the information-theoretical Araki-Lieb 
inequality [20], S < S ctl a 1 + S a2 a 2 , is maximized both by the two-mode squeezed 
vacuum and by a squeezed-coherent state, in agreement with the properties of 
one-mode squeezed states [18]. 
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According to eq.(19), the semiclassical entropy associated to a pure squeezed 
state \z) is thus non- vanishing, in spite of the fact that p = \z)(z\ = p 2 and S q (p) = 
0. This occurs because an effective reduction of the quantum density matrix is 
implicit in the semiclassical limit, as even a pure state is seen, classically, as a 
disordered one. However, the coarse graining associated to the semiclassical limit 
cannot give us explicit information about the physical origin of the decoherence 
mechanism, which can only be described in the context of some dynamic model 
for the quantum to classical transition (see also [21]). For a squeezed state such 
a model must account, in particular, for the fact that some degree of freedom 
of the system is made dynamically less relevant than others for what concerns 
macroscopic observations. 

Let us consider now the possibility of defining also a von Neumann entropy 
for the squeezing process, in terms of a properly reduced density operator. In the 
approach of Refs.[5-7] the reduction was performed by neglecting the off-diagonal 
elements of the density matrix in a superfluctuant representation. The use of 
such a representation is justified because the squeezed density operator can only 
be reconstructed, observationally, by measuring superfluctuant variables and their 
momenta. But is the reduction a formally justified procedure? 

The answer is certainly not for the pure squeezed vacuum state (9), whose 
density matrix in the superfluctuant representation 

p z (x, x', x, x') = (xx\z) (z\x'x') = ip z ( x j x )i>t( x ' ' i x ') = 




is perfectly symmetric in the (x, x') and (x, x') planes. The same answer holds 
for a pure squeezed number and squeezed-coherent state. On the contrary, the 
reduction would be naturally justified if the off-diagonal matrix elements would be 
suppressed, \p z ( x , x)\ » \p z (x, x' ^ x)\ (and the same for x), for any given value 
of the squeezing parameter. It is thus interesting to note that such requirement 
can be satisfied when the initial state is not pure, but is represented by a mixture, 
whose statistical weights tend to be equally distributed among all the states of the 
mixture. 

Suppose to start, in fact, from an initial configuration more general than the 
vacuum, represented by an ensemble of states which, like the vacuum, minimize 
however the quantum fluctuations Ax, Ax ( and their conjugate Ay, Ay). This 
corresponds to an initial density operator which can be expressed in terms of the 
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coherent states (15), (16) as 

Pi = J d 2 ad 2 aP(a, a)\aa){aa\ (23) 
where the statistical weights P > satisfy the normalization condition 

d 2 ad 2 aP(a 7 a) = 1 (24) 



The squeezed evolution leads then to a final density operator, pf = EpjEt, which 
in the superfluctuant basis is explicitly represented by 

Pf(x, x', x, x') = J d 2 ad 2 aP(a, a)ip za &(x, x)ip* a& (x', x) (25) 

Here ^^(i, x) is the squeezed-coherent wave function, representing in the (x,x) 
space the eigenfunctions of {a, a} with eigenvalues {a, a}, namely 



= -J -exp 



^zota{x,x) = (xx\zaa) = (xx\T,\aa) = 



-Umf - \0t) 2 + V^(ka + ~kx) - G -{x 2 + x 2 ) 



(26) 



(k,k are defined in eq.(17)). Note that, for pi ^ p 2 , there is a non-vanishing 
quantum entropy, S q = — J d 2 ad 2 aP{a 1 a)h\P{aa), naturally associated to the 
initial mixture. Such an entropy, however, is constant throughout the whole pro- 
cess of squeezing pi — > pf (which describes a unitary evolution), so that it cannot 
contribute to the overall entropy difference AS = Sf = S{. 

The choice of the vacuum as initial condition corresponds to the assumption 
P(a,a) = 5 2 (a)5 2 (a) in eq.(25). We shall consider here an initial distribution 
which, instead of being infinitely peaked upon the vacuum, tends to be spread 
rather uniformly over a wide range of states, so as to reflect our lack of knowledge 
about the initial conditions. We shall assume, however, that such a distribution is 
always centered around the vacuum, and we shall conveniently translate the (x, x) 
frame in such a way that (zaa\x\zaa) = = (zaa\x\zaa) (which implies 3?/c = 
= $lk, see eq.(26)). The simplest example satisfying the above requirements is 
a normalized square-like distribution, in which the probability of a state \aa) is 
constant and equal to 1/L for \a\, \a\ < L/2, and vanishing for \a\, \a\ > L/2, 
namely 

d 2 ad 2 aP(a 7 a) = da±da2daida2S(ai — a{)5{a.2 + 0^2) x 
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x -4 [0((*i + L/2) - 9{ ai - L/2)} [6(a 2 + L/2) - 6{a 2 - L/2)\ (27) 
L z 

Here 6 is the Heaviside step function, and the two delta distributions have been 
inserted to implement the conditions (x) = = (x); we have also rescaled the 
variables a, a so as to absorb the overall phase e~ ld which does not contribute to 
the integration over the whole complex plane. 

With this probability distribution the squeezed matrix (25) becomes 

p f {x,x,x,x)- ^ L ^ {x ,_ x) L ^_ £ /) e W 

and it is immediately evident that, for large enough L, the off-diagonal terms 
x 7^ x' , x 7^ x' are suppressed with respect to the diagonal ones. Such a suppression 
grows with L, namely with the number of states which are considered equiprobable 
in the initial distribution, i.e. with the growth of our "ignorance" about the initial 
configuration. For large L we may thus approximate the squeezed matrix with 
the diagonal elements only, and in particular we can formally define in this basis 
a (normalized) reduced density matrix, in the limit L — > oo (at finite a), as 

Pred(x, x', x, x') = lim {L 2 apf {x, x\ x, x) } = 

L^oo 

= - e -^ x2++ ^5(x - x')5(x - x') (29) 

TV 

This is exactly the representation in the (x,x) basis of the reduced operator 

Pred = J dxdi \M^)\^)^\ (30) 

previously used in Refs.[6,7]. 

Note that Trp re d = 1, but p 2 red ^ p re .d- This reduction leads then to a quan- 
tum entropy S q (p re d) which, unlike the entropy of the mixture S q (pf), depends on 
a, and in particular grows with the degree of squeezing. By taking the difference 
between final and initial entropy we find indeed the result of eq.(2), 

AS = (Trp red lnp red ) a=1 - (Trp red In p re d) a = 2r (31) 

The accuracy of this estimate in our context becomes larger and larger as the 
initial probability distribution becomes flatter, so as to share probabilities among 
a larger number of states. It is also important to note that if, for a given L, the 
reduced matrix represent a good approximation to the exact matrix (28) at large 
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squeezing (a >> 1), then the approximation is even better as we go back in time 
towards the initial configuration, where a = 1. 

We may thus conclude that the effective reduction of the squeezed density 
operator, introduced "ad hoc" in previous papers [6,7], can be formally justified 
through a limiting procedure based on an appropriate mixture of coherent states, 
assumed to represent the initial configuration. The corresponding entropy growth 
may thus physically interpreted as an amplification, due to the squeezed evolu- 
tion, of our ignorance about the initial conditions, in a context in which we have 
no motivation to single out some preferred initial state, and we cannot reconstruct 
the initial state through final measurements of superfluctuant variables only. Such 
an interpretation seems to be particularly appropriate in the context of cosmolog- 
ical perturbation theory, where the amplification of the background fluctuations is 
computed starting from an initial state which is left unspecified except for the re- 
quirement of minimizing the quantum fluctuations, and which can thus be properly 
represented as a general mixture of coherent states. 
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